Memory and aging effects in interacting sub-lOnm nanomagnets 
with large uniaxial anisotropy 



o 
o 

(N 



C/2 



X3 

o 



> 
o 
o 
m 

o 

On 
O 



X 



Kai-Cheng Zhang and Bang-Gui Liu 

Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China and 
Beijing National Laboratory for Condensed Matter Physics, Beijing 100190, China 

(Dated: June 6, 2009) 

Using a nonequilibrium Monte Carlo metliod suitable to nanomagnetism, we investigate represen- 
tative systems of interacting sub-lOnm grained nanomagnets with large uniaxial anisotropy. Various 
magnetization memory and aging effects are found in such systems. We explain these dynamical 
effects using the distributed relaxation times of the interacting nanomagnets due to their large 
anisotropy energies. 
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I. INTRODUCTION 



Nanomagnets attract huge interest because of their 
amazing properties and promising applications [ll, 0, H, 
01 S For well-separated nanomagnets (including 

single- molecule magnets), quantum tunnelling, inter- 
ference, and coherence can be observed at extra-low 
temperatures @, 0, H, H, and the magnetization be- 
haviours at some higher temperatures can be described 
by Neel-Brown law [ill \v\ . When inter-particle dis- 
tances become small enough, the dipolar-dipolar inter- 
action will modify the magnetization behavior leading to 
some super-spin-glasses behavior s [l3l [l^ , similar to con- 
ventional spin glasses [H, [l^ [13, [13] ■ Grained nanomag- 
nets with large uniaxial anisotropy are essential to mod- 
ern magnetic data storage. For typical CoCrPtB media, 
usual average grain sizes must be about lOnm to keep 
magnetic stability of 10 years at room temperature [l9l| . 
In the case of CoCrPt-oxide media for perpendicular 
recording, dominant inter-grain interactions are weak an- 
tiferromagnetic (AFM) couplings and average grain size 
can be 8nm or smaller for the same stability [l9j. Such 
average size can even be reduced down to 3nm or smaller 
when FePt in the Lip phase is used as data storage me- 
dia, because its magnetocrystalline anisotropy reaches to 
44 meV/nm^ [l9| . Such nanomagnets with large uniaxial 
anisotropy, especially when composing special systems, 
can yield various dynamical phenomena waiting for ex- 
ploration. 

Here we explore dynamical magnetic properties of rep- 
resentative systems of sub-lOnm grained nanomagnets 
with large uniaxial anisotropy. We use the giant spin 
approach for the nanomagnet because the magnetic in- 
teractions between electronic spins in it are strong. We 
assume that the magnetic anisotropy energies satisfy a 
Gaussian distribution to consider their fluctuations due 
to different shapes, sizes, and interfacial environments, 
and the inter-nanomagnet interactions, including mag- 
netic dipolar interaction, are AFM[i^. We use a giant- 
spin model and a dynamical spin Monte Carlo (DSMC) 
method [20l| to correctly simulate dynamical magnetiza- 
tion of the giant-spins of the component nanomagnets. 



Through systematical DSMC simulations, we find vari- 
ous field-cooling (FC) and zero-field-coohng (ZFC) mag- 
netization memory and aging effects in such systems. We 
explain these dynamical effects uniformly in terms of the 
continuously-distributed relaxation times of the compo- 
nent nanomagnets due to the various uniaxial anisotropy 
energies. More detailed results will be presented in the 
following. 



II. MODEL, METHOD AND PARAMETERS 

We consider a finite two-dimensional lattice of grained 
nanomagnets. Each nanomagnet actually includes many 
magnetic atoms, but the magnetic interactions between 
the magnetic atoms are much stronger than those among 
different nanomagnets. The inter-nanomagnet magnetic 
interactions are described by an antiferromagnetic cou- 
pling according to actual materials for modern magnetic 
data storage [l9l|. We take the giant spin approach and 
use one spin variable Si to describe the magnetic property 
of each nanomagnet. Because the spin value Si of such a 
nanomagnet is typically 10^ ^ 10'^, large anisotropy ener- 
gies of such nanomagnets are mainl y d ependent on their 
shapes and interfacial environments [ 1 9j . and thus can be 
reasonably assumed to satisfy a Gaussian distribution. 
This simplification keeps the main physics of these nano- 
magnet systems. Because S is large enough, we treat the 
spin operator Si as SiSi, where Si is a classical unit vector. 
Generally speaking. Si should vary from one nanomagnet 
to another, but for our systems all the {Si} are uniform 
enough to be let have the same average value S because 
the deviations from S make little differences. Therefore, 
our model can be described by the Hamiltonian, 



(1) 



where kui describes the uniaxial anisotropy energy satis- 
fying a Gaussian distribution with the average value 
and width 7 is defined as g^ofisS, and the field B is 
in the easy axis. Here Jij describes the inter-nanomagnet 
AFM interactions. For actual perpendicular media for 
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modern magnetic data storage, we assume that the easy 
axis is perpendicular to the plane of the N x N nano- 
magnet lattice. In such a setup, the magnetic dipolar 
interaction is reduced to an AFM intcr-nanomagnet in- 
teraction, which has been included in the Jy parameters 
in Eq. (1). 

Each of the spins has two meta-stable orientations 
along the easy axis and needs to overcome an en- 
ergy barrier to achieve a reversal. We use the DSMC 
method to simulate the spin dynamics of these nano- 
magnet svstems(20j. This method originates from the 
kinetic Monte Carlo method for simulating atomic dy- 
namics during epitaxial growth [21j]. Using 6i to describe 
the angle deviation of Si from the easy axis, we ex- 
press the energy increment of the i-th nanomagnet as 
AE, 
h,, = 



takes either 1 or -1[20]. As a result, an energy barrier 
AEi = {2kui + hiY /Akui must be overcome to achieve 
the reversal of the i-th spin. The rate for the spin re- 
versal obeys Arrhenius law Ri = Roexp{—AEi/kBT), 
where ks is Boltzmann constant and T is temperature. 

In our simulation, we use a typical value 1.0 x 10^/s 
for the characteristic frequency Rq. As for the anisotropy 
energy parameters, we reasonably assume that the aver- 
age value ku of {k^i} is 80.0 meV and the Gaussian width 
<7u is 44.7 meV. We change temperature T by a step 0.2 
K and set the sweeping rate to 1 K/s in most of the fol- 
lowing cases. Special cases will be explicitly described 
otherwise. We take = 40 and use periodic boundary 
condition. This lattice size is appropriate considering the 
actual situation in the media for the modern magnetic 
data storage. It almost does not matter which boundary 
condition is chosen with the lattice size. Further simu- 
lations with larger N are done for confirmation. Each 
of our data is obtained by calculating the average value 
over 500^^1000 independent simulation runs. 



III. MAIN SIMULATED RESULTS 

FC and ZFC magnetization curves. We simulate the 
FC curves by calculating the average magnetizations un- 
der 100 Oe at every T point when cooling from 110 to 10 
K. The ZFC curves are simulated by letting the system 
cool under zero field from 110 to 10 K and then calculat- 
ing the average magnetizations under 100 Oe at every T 
point when warming from 10 to 110 K. Our simulation re- 
alizes what happens in measuring experimental FC and 
ZFC magnetization curves. Our simulated results are 
presented in Figure 1. For zero interaction or J = 0, 
the FC magnetization increases monotonously with de- 
creasing T, and the ZFC magnetization increases first 
with increasing T, reaches its maximum at the block- 
ing temperature Tg = 55 K, and then decreases. The 
FC and ZFC magnetization curves follow the Curie law 
above the temperature Tm = 70 K. The difference be- 
tween the FC and ZFC magnetizations diminishes above 



Tm- Such behaviors are key features of nanomagnets, 
showing up in su per -spin glass systems [3] and interact- 
ing nanoparticles[22. [23j . When J is larger, the magneti- 
zation in the FC and ZFC curves becomes substantially 
smaller, as shown in Figure 1. When J becomes further 
larger, a minimal magnetization can be seen below Tb- It 
is at 30 K for 2.0 meV. When J is less than 0, the mag- 
netization in such curves becomes larger but the curve 
shape remains nearly the same. 



kui sin 9i — hi{cos9i — 1) to leading order, where fo 
{jB — '^jJijSj)si and the reduced variable Si 
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FIG. 1; Simulated FC (black) and ZFC (gray) magnetization 
curves under a field 100 Oe for J = 0.0, 1.0, and 2.0 meV. 
These curves are calculated in the same way as corresponding 
experimental magnetization curves are measured. 



Simulated ZFC and FC memory effects over time. Cool 
the system under zero field to a low temperature, 30 K, 
and then apply a field 500 Oe and let the system relax 
from i = 0. At ti — 814 s, we change T to 10 K and re- 
verse the field, and at t2 — 4124s, we recover the original 
temperature and field. We calculate the average magne- 
tization when the system relaxes with time from t = Q 
to t — 15000s, and present the results in Figure 2a. The 
magnetization at t2 is equivalent to that at ti, that is, 
the system keeps the memory at ti although it undergoes 
relaxation from ti to t2 under the different temperature 
and field. This is a ZFC memory effect over time. In 
addition, we cool the system under a field 500 Oe to 30 
K, and then remove the field and let the system relax 
from t = 0, but we change the temperature to 10 K and 
the field to -500 Oe between t[ = 3416s and t'^ = 6730s. 
Meanwhile we calculate the magnetization from t — to 
t = 20000s. The result, shown in Figure 2b, shows an 
FC memory effect over time, mf^ = . Similar effect 
was observed in super-spin glass systems and interacting 
nanoparticlesp^ 123. [23 1 . Our further calculations show 
that a small inter-spin interaction, either AFM or FM, 
does not affect the effects substantially. 

Simulated FC memory effect over temperature. Cool 
the system under a field 100 Oe from 110 to 10 K, but 



3 



0.56 (a) _^j=omeV 
— ^ J=1.0 meV 



c 0.42 
o 

to 

N 

1 0.28 - 



0.14- 



10 K 

-500 Oe 




30 K 
500 Oe 



30 K 
500 Oe 



0.00 L 





(b) 


-J=0 meV 






-J=1.0 meV 




30 K 






OOe 








30 K 




OOe 


-500 Oe 

L 





4 8 12 16 5 10 15 20 



Time (10' s) 



Time (10' s) 



FIG. 2: Simulated ZFC (a) and FC (b) memory effects over 
time with J — and 1.0 meV. (a) The temperature is 30 K 
and the field 500 Oe at all the time points except between 
ti = 814 s and t2 = 4129 s (gray); in the special time period, 
the temperature and field are changed to 10 K and -500 Oe. 
(b) The temperature is 30 K and the field Oe at all the time 
points except between t'l — 3416 s and t2 = 6730 s (gray); in 
the special time period, the temperature and field are changed 
to 10 K and -500 Oe. 
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FIG. 4: Simulated ZFC memory effects over temperature for 
J — 1.0, 1.5, and 2.0 meV. The references (black lines) are 
standard ZFC magnetization curves with 100 Oe, and the 
memory curves (gray dash lines) are made by keeping the 
system relaxing for additional 100 s at 40 K during the cool- 
ing. The heating is always homogeneous in time and field. 
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FIG. 3; Simulated FC memory effects over temperature for 
J = 0.0, 1.0, and 2.0 meV. During cooling (black), the tem- 
perature is changed homogeneously and the field is 100 Oe at 
all the temperature points except 40 K and 20 K. At these two 
special temperature points, the system is kept at the temper- 
ature for additional 40 s and 60 s, respectively, and the field 
is kept constantly the same -100 Oe meanwhile. The heating 
(gray dash) is homogeneous in time under the same field. 



at two temperature points, 40 K and 20 K, the field is 
reversed and then recovers after 40 s and 60 s, respec- 
tively. The resulting FC magnetization curves for J = 0, 
0.5, and 1 meV are shown in Figure 3. The step is clearly 
seen at both 40 K and 20 K for each of the three curves. 



The abrupt magnetization change is because the field is 
kept as the reversed value, -100 Oe, for additional 40 s 
or 60 s. When the cooling ends at 10 K, we warm the 
system under 100 Oe from 10 to 110 K. The correspond- 
ing three warming magnetization curves are presented 
in Figure 3. Clear magnetization drop is observed at 
both 40 K and 20 K, although the warming is homo- 
geneous in time and field. This means that the system 
during the warming has the memory of the abrupt mag- 
netization changes at the same temperatures during the 
cooling. Similar effects have been observed in interact- 
ing nanoscale systems [T^. [2^ [23j . Our results show that 
a small inter-spin interaction changes the FC memory 
effect only a little. 

Simulated ZFC memory effect over temperature. The 
above effects are essentially independent of inter-spin in- 
teractions J. Such interactions weaken, even can break 
these memory effects, but appropriate J values can lead 
to another memory effect shown in Figure 4. First, cool 
the system homogeneously under zero field, calculate 
standard ZFC magnetization curves when warming the 
system under 100 Oe, and take them as references; then, 
cool the system under zero field but keep it for additional 
100s when T = 40 K, calculate corresponding magneti- 
zation curves when warming the system under 100 Oe. 
Although the warming is homogeneous in time, the mag- 
netization near 40 K is smaller, which implies that the 
system has the memory of the waiting at 40 K during 
the cooling. These are consistent with observations in 
the cases of other interacting super-spin systems [l^, H^- 
It is clear that nonzero J values are helpful for the mem- 
ory effect. 
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FIG. 5: Simulated aging eflect. The magnetization m (a) and 
rate S (b) as functions of time t at 15 K under 100 Oe for 
three waiting time periods t^,: 100s (solid), 1000s (dash), and 
10000s (dot). The inset in (a) shows the same m results in 
the linear scale of t. 



Simulated aging effect. Cool the system under zero 
field to 15 K (below Tb) and keep it unchanged for a pe- 
riod i^o, and then let it relax with time t under 100 Oe. 
During the relaxation, we calculate the magnetization 
m and the rate S = dm/d\ogt as functions of t. Our 
simulated results in logarithmic scale with t^ — 100s, 
tw = 1000s, and t^ = 10000s are presented in Figure 5. 
The inset in the upper panel shows the magnetization in 
linear scale of t. It is clear that the longer the waiting 
time, the smaller the magnetization, and the relaxation 
rate has a maximum sX t k, tw This is consistent with 
well-known aging effect observed in spin-glass and inter- 
acting super-spin systems [H, [H, [13, ll^ • 



IV. DYNAMICAL MECHANISM AND 
COMPARISON WITH EXPERIMENT 



of nanomagnets. The key point is that for a given time 
period and at a given temperature, most effectively re- 
laxed are the nanomagnets with appropriate relaxation 
time T because those with smaller r have been sufficiently 
relaxed and those with larger r remain frozen. For the 
memory effects over time shown in Figure 2, the temper- 
ature between the two time points, ti and t2 (or t'l and 
^2), is substantially lower, and nothing important hap- 
pens meanwhile, although the field is changed, because 
the nanomagnets with much smaller t have already been 
thoroughly relaxed, and hence the magnetization is re- 
covered after the field and temperature are restored. As 
for the memory effect shown in Figure 3, the additional 
relaxation under the reversed field at both 40 K and 20 
K during the cooling makes corresponding nanomagnets 
orient in the reversed direction and most of them, re- 
maining frozen at the lower temperatures in the further 
cooling, cause the magnetization dips at both 20 K and 
40 K in the warming curve, the FC memory effect over 
temperature. For the ZFC memory effect shown in Fig- 
ure 4, the additional 100 s relaxation at 40 K during the 
cooling, in the presence of the AFM inter-nanomagnet 
coupling, makes the nanomagnets with corresponding 
anisotropy energy over-relaxed so that the magnetiza- 
tion curve (warming from 10 K) has a small dip at the 
same 40 K with respect to the reference curve. As for 
the aging effect shown in Figure 5, the key point is that 
the system is still relaxing at the base temperature 15 K 
during the additional waiting period t^. For larger t^, 
the system becomes more relaxed and the magnetization 
is a little smaller until a time t Kt^. However, the effect 
of the additional relaxation of i^, will finally disappear 
for t ^ tu,- Therefore, there is a peak for the rate of 
increasing the magnetization at a time t k, tw 

Compared with experiments and other theories. Sim- 
ilar FC and ZFC magnetization memory effects were 
observed in both spin (including super-spin) glass sys- 
tems and isolated nanomagnets, and similar ZFC mem- 
ory and aging effects in spin glass systems [l^ [13, 
El [S EI 13, 113. Usually, memory effects in iso- 
lated nanomagnets were attributed to broad distribu- 
tion of relaxation times for different nanomagnets. The 
memory and aging effects in spin glass systems can 
be understood in terms of droplet model and hierar- 
chy model [ij, II6I, im. Ill- 

For our systems the mag- 
netic anisotropy is uniaxial and the inter-spin interaction 
is uniformly AFM, and therefore our model, including 
inter-spin interaction, is dis tinguished from other ones 
for those systems 
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Mechanism for the memory and aging effects. Our sys- 
tems consist of nanomagnets with large uniaxial magnetic 
anisotropy which follows Gaussian distribution. The dis- 
tribution of anisotropy energy causes various relaxation 
times through relation = tq exp(A£'i/fcBT). The 
memory and aging effects can be explained in a unified 
way by the continuously-distributed relaxation times {ri} 



V. CONCLUSION 

Using the dynamical spin Monte Carlo method, we in- 
vestigate various nonequilibrium dynamical magnetiza- 
tion behaviors of the representative systems (finite two- 
dimensional lattices) composed of antiferromagnetically- 
coupled sub-lOnm grained nanomagnets with large uni- 
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axial anisotropy. In each of such systems of the nanomag- 
nets with the same easy axis, we find both ZFC and FC 
magnetization memory effects over both time and tem- 
perature and the magnetization aging effect, which partly 
appeared otherwise in different spin (or super-spin) glass 
systems and interacting nanoparticles. We explain these 
interesting dynamical effects in a unified way in terms 
of the continuously-distributed relaxation times of the 
interacting nanomagnets due to their large anisotropy 
energies in the actual media. 
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